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ABSTRACT 

The  purpose  of  the  study  was  to  further  understanding  of  unsteady  transport  of  coarse  bedload.  A 
particular  focus  was  the  interaction  of  the  load  with  the  river  bed  material  and  the  effects  of  transport  on 
stream  hydraulics.  It  was  anticipated  that  new  methodology  would  be  developed  and  that  data  sets 
suitable  for  further  modelling  would  result.  The  nature  of  coarse  gravel  transport  through  a  short  reach 
of  a  mountain  river  was  examined  using  novel  methodology;  primarily,  magnetic  particle  detection, 
acoustic  bedload  detection,  electromagnetic  current  meters  and  fine-resolution  bed  elevation  surveys. 

For  a  wide  range  of  flows,  up  to  and  including  bankfull,  the  surface  of  an  alternate  bar  and  the 
neighbouring  talweg  channel  remained  essentially  stable.  At  a  critical  in-bank  discharge,  bed  sediment 
was  transported  from  sources  upstream  through  the  study  teach  as  a  discrete  slug  of  bedload.  The 
interaction  of  this  mobile  ’sheet’  with  the  coarse  bed  sediments  led  to  selective  deposition  of  fine 
bedload  into  the  interstices  of  the  static  armour.  This  mechanism  caused  physical  smoothing  of  the  bed 
and  concomitant  reduction  in  hydraulic  roughness  during  the  passage  of  the  sheet.  Consequent  energy 
slope  adjustment  caused  water  depths  to  drop  and  velocities  to  increase,  leading  to  local  readjustments 
in  shear  stress  and  streampower.  As  the  tail  of  the  sheet  passed,  the  fine  particles  were  reentrained  from 
the  armour,  the  bed  roughened  and  hydraulic  parameters  adjusted  again. 

It  is  concluded  that  although  exceeding  an  hydraulic  threshold  is  required  to  initialize  bedload  transport, 
there  are  important  mutual  feedback  mechanisms  whereby  the  interactions  of  the  mobile  sediment  and 
the  bed,  lead  to  further  adjustment  in  the  hydraulic  climate  which  then  further  affect  bedload  transport. 
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INTRODUCTION 

The  process  of  bedload  transport  in  gnvel-bed  streams  is  of  practical  interest  to  engineers  and  attracts 
considerable  theoretical  and  empirical  attention.  One  important  focus  has  been  the  nature  of  particle 
interactions  with  the  bed  at  the  temporal-scale  of  turbulence  measurements.  At  a  larger  scale,  equivalent 
to  the  passage  of  an  hydrognph,  a  second  focus  has  been  the  hydraulic  conditions  necessary  for  initial 
transport  and  the  relationship  between  flow  and  the  total  quantity  of  material  in  motion  close  to  the  bed. 
Recently  limited  field  and  laboratory  data  have  indicated  that  bedload  motion  is  distinctly  unsteady  even 
during  hydrognphs  characterized  by  competent  flow  and  is  intimately  related  to  turbulent  flow  structure 
at  fine  spatial  and  temporal  resolution.  Bedload  may  move  as  discrete  pulses  resulting  in  intense  periods 
of  transport  and  vice  versa.  Concomitant  changes  in  bed  elevation,  bed  grain-size,  hydraulic  roughness 
and  near-bed  velocity  structure  are  expected.  Until  recently  there  were  no  comprehensive  data  sets  that 
could  be  used  to  demonstrate  the  relationship  of  bedload  pulses,  channel  boundary  and  hydraulic 
interactions  in  natural  rivers.  Data  available  relate  almost  exclusively  to  sand  and  fine  gravel.  Distinctly 
lacking  are  data  pertinent  to  the  motion  of  coarse  gravel  incl  -ding  pebbles  and  cobbles  up  to  200mm  in 
size.  The  purpose  of  this  study  was  to  develop  techniques  which  would  provide  comprehensive  field 
data  related  to  the  passage  of  sheets  of  coarse  bedload.  At  the  same  time  it  was  anticipated  that  daU  sets 
would  result  which  would  be  suitable  for  mathematical-modelling  purposes. 

Technical  Objective 

Rather  than  investigate  initial  motion,  particle  transport  or  bed  level  fluctuations  in  isolation  in  a  natural 
stream,  the  present  objective  was  to  monitor  the  bedload  phenomena  and  associated  processes 
intensively  over  the  full-width  of  a  defined  reach  where  the  inputs  and  outputs  of  sediment  were  known. 
Linkages  between  hydraulics,  initial  motion  criteria,  transport  mechanisms  and  bed-lcvcl  fluctuations 
were  to  be  established.  A  better  understanding  of  the  complex  processes  of  sediment  transport  and  river 
channel  adjustment  should  ensue.  Techniques  developed  and  data  acquired  should  have  immediate 
relevance  to  defining  theoretical  bedload  transport  relationships  and  models  of  sediment  sorting  which 
have  practical  application  to  hydraulic  engineering. 

Statement  of  Programme 

The  methods  are  fully  explained  below  but  a  brief  synopsis  at  this  point  is  useful.  The  program  was 
conducted  in  the  gravel-bed  Squaw  Creek  in  Montana,  USA,  near  the  city  of  Bozeman  between  1990 
and  1992.  The  river  is  ideal  for  study  as  it  has  a  snow-melt  regime  in  April  and  May  each  year.  Steady 
changes  in  stage  occur  each  day  in  response  to  diurnal  changes  in  snowpack  melt  rates  such  that 
competent  flows  and  bedload  transport  often  occur  at  the  study  site  in  the  evening.  Individual  naturally- 
magnetic  particles  and  kinematic  waves  of  particles  were  recorded  passing  into  the  experimental  reach 
and  exiting  from  the  reach  using  two  electromagnetic  coils  placed  across  the  river  bed.  Within  the  reach 
a  transect  was  established  across  a  section  which  included  the  main  talweg  anti  an  adjacent  lateral  bar. 
Detailed  bed  level  fluctuations  across  this  section  were  recorded  and  detailed  velocity  profiles  taken, 
from  which  shear  stress  distributions  across  the  bed  were  derived.  Bedload  samples  obtained  using  a  3* 
He  I  ley -Smith  and  a  large  net-sampler  were  used  to  calibrate  the  magnetic  coil  record.  These  methods 
had  been  developed  at  Squaw  Creek  prior  to  the  present  investigation  and  had  yielded  small  amounts  of 
useful  data  as  techniques  were  refined.  The  three  years  of  the  present  study  were  unusually  dry  and  only 
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one  season  resulted  in  bedload  motion.  Despite  this  disappointment,  equipment  was  install'  d  and 
improved  in  the  fust  year,  four  flood  waves  were  monitored  in  the  second  year  and  in  the  final  year 
turbulence  data  were  obtained. 

Historical  Background 

An  important  focal  point  to  river  research  in  recent  years  has  been  the  nature  of  transport  in  gravel-bed 
rivers  over  time-scales  comparable  to  the  passage  of  an  bydrograph  and  at  length-scales  equivalent  to 
the  channel-width.  Answers  to  a  variety  of  questions  are  currently  being  sought.  For  example,  what  are 
the  hydraulic  conditions  necessary  to  initiate  transport  of  natural  bed  mixtures  and  what  is  the  nature  of 
transport  as  an  event  begins  and  ends?  Even  during  steady  hydraulic  conditions,  it  has  long  been  known 
that  gravel  particles  may  move  in  groups,  forming  bedload  pulses1,  so  that  the  bedload  transport  rate  is 
unsteady2  and  also  may  vary  spatially  across  a  section.  As  pulsing  is  associated  with  scour  and  fill3,  it 
follows  that  the  passage  of  a  large-scale  pulse  should  be  manifest  as  a  change  in  the  width  or  depth  of 
the  stored  volume  of  bed  sediment;  inducing  the  development  of  dynamic  bed-waves2.  In  addition 
segregation  of  particles  by  size  within  pulses  may  occur  with  coarse  gravel  at  the  leading  edge  and  finer 
material  in  the  tail4. 

The  linkages  between  particle  motion,  bed  elevation  and  hydraulics  are  most  readily  addressed  in  the 
laboratory  flume2,3'6,  but  it  is  evident  that  success  in  obtaining  relevant  data  in  the  field  is  primarily 
constrained  by  technical  difficulties.  Most  notable  is  the  need  to  develop  continuous  automatic  detection 
systems  which  will  provide  adequate  temporal  and  spatial  resolution. 

In  the  past,  a  variety  of  devices  and  techniques  have  been  developed  to  obtain  suitable  field  data. 
Painted,  tagged  and  radioactively-  or  magnetically-labelled  clasts  have  been  used  as  tracers1,7 A'UQ,U. 
Mechanical  samplers,  such  as  the  Helley-Smith  sampler,  have  also  been  used  to  measure  coarse  bedload 
transport12,13.  Large-scale  pit  traps  (of  varying  complexity)  and  vortex  samplers  have  been  installed  in 
stream  bottoms14,13.  Instruments  have  been  developed  to  monitor  the  motion  of  individual  particles  on 
the  stream  bottom  as  well.  For  example,  in  1973  Tyoniuk  and  Wamock16  used  acoustic  techniques  to 
relate  the  sound  of  motion  of  individual  particles  on  the  stream  bottom  during  transport  to  velocity  and 
so  isolate  the  locus  of  threads  of  high  transport  intensity. 

The  above  methods  have  all  suffered  from  various  problems,  including  logistic  and  technical  constraints 
as  well  as  sampling  efficiency  and  avenging  problems.  The  list  of  techniques  is  not  exhaustive,  but  does 
serve  to  point  out  the  diversity  of  approaches  used  to  investigate  bedload  transport  in  gravel-bed 
streams.  The  proliferation  of  techniques  is  related,  in  part,  to  the  difficulties  associated  with  monitoring 
transport17,11.  Floodwaters  can  be  turbid  and  damaging  to  sensitive  equipment.  Adequate  sampling  is 
often  dangerous  to  the  researcher  because  it  must  be  done  during  peak  stage.  Natural  events  which  must 
be  sampled  occur  at  times  which  may  be  unexpected  and  may  occur  at  night.  These  problems  make  it 
essential  that  continuous  robust  automatic  monitoring  be  developed  with  the  time  and  spatial  resolution 
necessary  to  advance  bedload  science. 

In  recent  years,  potentially  suitable  detector  systems  have  been  installed  in  stream  beds  to  examine  the 
temporal  variation  in  (nab. rally-magnetic)  bedload  transport  setni-continuously11,19,  whilst  real-time 
acoustic  tracing  of  individual  particles  is  under  development.  The  present  generation  of  magnetic 
detector  systems  can  only  be  used  in  those  streams  where  the  local  geology  provides  a  ready-source  of 
naturally  magnetic  stones.  However,  in  those  situations  it  shows  great  potential  for  continuous  detailed 
monitoring  at  a  suitable  time  and  spatial  resolution. 
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Most  alternative  approaches  avenge  measurements  or  integnte  over  a  discrete  time  period  because  of 
mechanical  or  other  practical  constraints  .  Such  avenging  can  obscure  the  detail  of  transport 
patterns30,21'32.  For  example  if  10  minute  sediment  transport  avenges  are  used,  dynamic  features  with 
durations  of  less  than  10  minutes  are  lost22.  A  device  was  introduced  by  Ergenzinger  and  Custer19,  which 
potentially  reduces  such  problems  and  facilitates  the  investigation  of  detailed  transport  phenomena.  The 
instrument  is  described  below  and  is  capable  of  counting  the  passage  of  individual  naturally  magnetic 
pebbles  and  cobbles.  Consequently  the  device  produces  a  real  time  record  of  sediment  transport  with 
minimal  count-averaging.  Such  data,  when  combined  with  hydraulic  data,  make  it  possible  to  study  the 
nature  of  transport  as  coarse  particle  motion  begins  and  ends  and  to  examine  unsteady  sediment 
transport. 


Study  Site 


Squaw  Creek  is  situated  in  the  Northern  Rocky  Mountains  of  the  USA  and  is  a  tributary  of  the  Gallatin 
River  which  together  with  The  Madison  and  Jefferson  Rivers  joins  to  form  the  Missouri  approximately 
65km  downstream  from  the  Squaw  Creek  confluence  (Fig.  1).  Squaw  Creek  drains  an  area  of  106km2  of 
andesitic  volcanic  terrain  which  comprises  the  Gallatin  Range,  Gallatin  County,  Montana.  Originating 
below  Hyalite  Peak,  Squaw  Creek  flows  for  approximately  22km  through  a  forested  catchment  before 
joining  the  Gallatin.  The  climate  of  this  area  is  characterised  by  long,  cold  winters  and  short,  cool 
summers.  The  avenge  annual  precipitation  of  813mm  consists  of  a  snowfall  equivalent  of  457mm  and 
356mm  as  rain.  The  snowpack  in  the  mountains  reaches  its  maximum  water  content  in  April  or  May 
and  melting  is  usually  completed  by  June  or  July*.  Rapid  melting  has  been  observed  in  late  May  and 
early  June  in  recent  years  (Fig.  2)  and  the  river  generally  reaches  its  peak  annual  discharge  at  this  time. 
Peak  discharges  are  caused  typically  by  snowmelt  or  more  occasionally,  and  less  predictably,  as  a  result 
of  convective  rainstorms.  As  the  lag  time  for  snowmelt  is  about  8  hours,  the  peak  diurnal  stage  is  not 
reached  until  circa  2000  hr  when  discharge  is  typically  <6  m5  s  ’.  Lag  time  for  rain-induced  hydiugraphs 
is  less;  typically  1.5  hours,  commonly  giving  peak  flows  in  the  early  afternoon.  Bankful]  discharge  is 
about  7  mJ  s'1  with  a  fifty-year  flood  equal  to  20  m1  s'1  (Fig.  3).  Further  details  of  the  hydrological 
conditions  are  given  elsewhere34,25,2*27. 

The  bedrock  underlying  the  drainage  basin  of  Squaw  Creek  (Fig.  4)  is  rich  in  magnetic  minerals, 
composed  of  approximately  25%  Archaen  quartzo-feldspathic  gneiss,  20%  Palaeozoic  and  Mesozoic 
sediments  and  50%  Eocene  andesitic  and  basaltic  lava  flows,  mudflows  and  igneous  intrusions3*.  Further 
details  of  the  drainage  basin  characteristics  are  given  by  Ergenzinger  and  Custer24. 

The  30m  long  study  reach  is  situated  100m  upstream  of  the  confluence  with  the  Gallatin  river,  and 
consists  of  an  alternate  bar  on  the  left  bank  with  the  main  talweg  channel  against  the  right  bank. 
Immediately  upstream,  but  outwith  the  study  reach,  an  alternate  bar  exists  on  the  right  bank.  During 
very  high  flows  (>7mJ  s'1)  a  secondary  channel  is  inundated  on  the  right  bank  (Fig.  5).  This  channel  was 
rarely  active  during  the  study  period  with  negligible  sediment  transport  and  was  not  included  in  the 
present  study.  The  bed  slope  of  the  river  locally  varies  between  0.02  and  0.03,  whilst  the  channel  width 
varies  between  8m  and  20m.  Through  the  study  reach  the  bed-width  is  constant  (10m).  The  channel  has 
a  width  to  depth  ratio  of  about  6.7  at  bankful]. 

Median  grain  sizes  (DM)  of  the  bulk  bed  material  range  from  19mm  on  the  bar  to  140mm  in  the 
channel,  where  more  than  16%  of  the  particles  are  larger  than  200mm.  According  to  analyses  by  Bunte39 
the  majority  of  particles  are  ellipsoidal,  falling  in  the  central  bladed-class  of  Sneed  and  Folk30.  The  bed 
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sediments  (Table  1)  are  generally  well  armoured  by  material  >DM  of  bulk  samples;  that  is  90mm  on  the 
hart  op  and  200mm  in  the  base  of  the  channel.  This  armour  is  rarely  disrupted  although  flows  up  to  7  m3 
s  '  lead  to  slight  disturbance31  Bedload  up  to  180mm  has  been  recorded  for  Dows  up  to  bankfull31, 
indicating  that  the  larger  elements  in  the  bed  sediment  remain  as  a  stationary  lag.  These  large  particles, 
often  stand  proud  of  the  bed  level,  are  widely  scattered  over  the  bed,  and  may  provide  a  second  scale  of 
flow  resistance  in  addition  to  the  grain  resistance  imparted  by  the  periodically  mobile  finer  bar-top 
sediments. 

The  study  reach  exhibits  a  number  of  attributes  which  contribute  to  it’s  suitability  for  research  on  coarse 
bedload  transport  dynamics.  Most  importantly  Squaw  Creek  is  an  active  gravel-bed  stream  whose  banks 
are  stable  but  whose  bed  can  adjust  to  naturally  occurring  changes  in  sediment  supply  and  stream 
discharge.  The  diurnal  peak  discharges  in  Spring  are  responsible  for  almost  all  annual  bedload  transport. 
The  timing  of  these  hydrographs  are  more  or  less  predictable  and  commonly  have  similar  durations  and 
peak  values.  This  allows  the  collection  of  replicate  flood-event  data  during  a  defined  sampling  season. 

In  the  early  stages  of  the  development  of  the  magnetic  pebble  detector24,  the  Squaw  Creek  catchment 
was  identified  as  having  a  suitable  proportion  of  magnetic  bed  material  larger  than  30mm29  required  for 
detection.  Prior  to  the  present  study,  a  research  project  was  conducted  intermittently  between  1981  and 
1988  by  the  Free  University  of  Berlin  at  the  same  site.  These  studies  have  provided  background 
information  on  the  fluvial  geo  morphology  of  Squaw  Creek. 

2.  INSTRUMENTATION  AND  SAMPLING  METHODS 

In  summary,  the  instrumentation  consists  of  two  magnetic  detector  logs  (which  are  embedded  flush  with 
the  river  bed)  at  each  end  of  the  study  reach  and  a  measuring  platform  which  was  used  to  collect 
hydraulic  data  and  information  on  bed  level  variations.  Five  stage  recorders  have  been  installed  (Fig.  S) 
and  a  hut  erected  to  protect  instrumentation  and  provide  shelter  during  fieldwork.  Vehicular  access  is 
easy  and,  in  addition,  power,  flood-lighting  and  a  telephone  were  also  installed. 

a)  Detection  of  bedload  transport  using  the  magnetic  pebble  detector  system. 

The  particle  movement  detectors  (Fig.  6)  installed  at  Squaw  Creek  record  the  passage  of  individual 
magnetic  particles  both  spatially  and  temporally.  The  method  of  detection  of  particle  motion  is  based  on 
the  electromagnetic  Faraday  principle.  If  a  naturally  magnetic  or  artificially  magnetized  particle  is 
moved  over  an  electric  coil  the  motion  induces  a  voltage  fluctuation  in  the  coil.  This  voltage  is  produced 
by  induction,  and  consequently  no  electrical  power  is  supplied  to  the  coil.  However  amplification  of  the 
signal  is  needed  so  that  it  may  be  recorded.  During  the  development  of  the  prototype  system2332, 
cobbles  either  were  implanted  with  magnets  and  then  placed  on  the  stream  bed  or  naturally  magnetic 
clasts  were  used.  Owing  to  the  reduced  magnetism  of  the  natural  tracer  compared  to  that  of  the  artificial 
tracer  it  was  necessary  to  place  the  detector  coils  flush  with  the  stream  bed  surface  rather  than  bury 
them.  The  use  of  the  naturally  magnetic  tracer  was  superior  to  the  artificially  magnetized  clast  not  only 
from  a  practical  angle  but  their  use  also  eliminates  the  logistic  problems  associated  with  seeding 
artificial  tracers  on,  or  in,  the  stream  bed. 

At  Squaw  Creek  a  magnetic  particle  motion  detector  had  already  been  installed  at  the  downstream  end 
of  the  study  reach  by  the  Berlin  team.  During  the  present  study  period  an  additional  detector  was 
positioned  at  the  upstream  limit  of  the  reach  (Fig.  4).  Each  detector  system  has  the  capacity  to  detect 
motion  at  discrete  l.Sm  sections  transverse  to  the  flow.  This  configuration  was  devised  so  that  input  and 


Page  14  of  46  Pages 


DAJA45-90-C-0006 
FINAL  REPORT 


output  of  particle  numbers  to  and  from  the  reach  could  be  determined.  The  detectors  record  the  motion 
of  particles  in  real-time  and  their  signals  presently  are  counted  and  logged  at  10  minute  intervals.  The 
early  system  used  by  the  Berlin  team  often  was  unreliable,  with  false  signals  induced  by  static  electricity 
in  the  atmosphere,  for  example,  or  individual  signals  could  not  be  isolated  at  high  transport  rates,  when 
noise  to  signal  ratios  were  high*.  Logging  consisted  of  a  chart  record  which  had  to  be  manually 
analyzed.  Technical  information  on  this  prototype  detector  is  given  elsewhere33.  However,  during  this 
project,  considerable  technical  improvements  were  made  to  the  amplification/filter  system.  The 
electronics  consisted  of  four  main  parts,  a  double  auto  zero  amplifier  and  analogue  multiplexer,  a  data 
acquisition  controller  with  digital  converter  and  some  auxiliary  electronics  with  a  micro-controller.  This 
meant  that  the  signal  to  noise  ratio  was  improved  considerably.  Data  were  recorded  directly  on  a  PC 
with  visual  output  for  monitoring  the  signals. 

b)  Detection  of  bedload  transport  using  other  samplers 

Previous  research  at  the  site  had  made  use  of  hand-held  3*  or  6”  Helley-Smith  bedload  samplers  but 
these  had  proved  to  be  inadequate  in  high  Dow  and  incapable  of  trapping  the  larger  fraction  of  bedload 
in  transport28.  Latterly,  a  direct  method  of  bed  load  measurement  was  used  at  the  downstream  limit  of  the 
study  reach  by  trapping  all  gravel  greater  than  11mm  in  a  specially  designed  net  supported  on  a  frame31. 
This  large  bedload  sampler  was  deployed  below  a  sill  immediately  downstream  of  the  lower  detector 
(Fig.  4).  Although  close  in  proximity,  the  two  systems  did  not  interfere  with  each  other,  and  direct 
sampling  allowed  validation  and  calibration  of  the  detector  system29-31.  The  bedload  net  sampler  was 
deployed  for  variable  lengths  of  time  dependant  upon  the  prevailing  transport  conditions.  In  addition,  in 
1992  a  vortex  sampler11  was  constructed  to  obtain  bedload  bulk  samples  every  5  minutes,  but  no 
significant  transport  occurred  during  the  sampling  campaign. 

A  hydrophone34  was  calibrated  to  record  the  sediment  generated  noise  (SGN)  arising  from  inter-particle 
collisions  during  transport  of  bedload  material.  Unfortunately,  and  unusually,  no  bedload  transport 
occurred  in  the  1992  season  owing  to  inadequate  snow-pack.  However,  the  system  was  calibrated 
through  artificial  production  of  SGN.  This  passive  acoustic  method  provides  an  opportunity  to  obtain 
continuous  high  temporal-resolution  records  of  mass  transport  rates  with  negligible  interference  with  the 
state  of  the  bed  or  the  flow. 

The  idea  of  utilizing  noise  generated  by  bedload  to  quantify  transport  appears  to  have  originated  with 
MuhlhofTcr13.  Previous  attempts  to  apply  SGN  to  monitor  bedload  transport16-35,3<U7'H,0<Q'<l  u<3'44'45  have 
met  with  only  limited  success.  Using  a  combined  theoretical  and  experimental  approach  Thome46 
showed  the  SGN  spectrum  to  be  composed  of  band-limited  noise,  with  a  frequency  range  inversely 
proportional  to  the  mass  of  mobile  material.  The  origins  of  SGN  have  subsequently  been  explained 
theoretically  through  application  of  rigid  body  acoustic  radiation  theory47. 

Improved  understanding  of  marine  bedload  transport  processes  has  resulted  through  combined 
measurements  of  flow  turbulence  and  SGN.  In  recent  studies  conducted  in  the  west  Solent,  UK,  quasi- 
fluid  bursting  events  analogous  to  sweeps  were  found  to  be  correlated  strongly  with  intermittent  bedload 
transport  events34,48,49  and  improved  description  of  mobilization  threshold  and  transport  rates  have 
resulted.  More  recently,  attempts  have  been  made  to  monitor  mobilization  and  transport  of  gravel  under 
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combined  wave  and  current  conditions  (Hardisty,  personal  communication)  and  statistical  simulation 
studies  of  acoustic  bedioad  data  have  been  reported  by  Williams  and  Tawn30. 

Given  the  success  in  utilizing  SGN  and  EMCM  turbulence  measurements  in  past  marine  studies,  in  the 
present  study,  an  attempt  was  made  to  apply  similar  techniques.  This  is  the  first  attempt  to  deploy  the 
instrumentation  in  shallow,  fast-flowing  rivers. 

c)  Hydraulic  measurements 

i)  Water  depth:  Five  temporary  float-operated  water-level  recorders  were  deployed  in  stilling  wells 
positioned  in  the  stream  against  the  stream  banks  (Fig.  5).  The  chart  records  from  these  proved  useful 
for  routine  assessment  of  the  direction  of  changes  in  discharge  during  individual  hydrographs.  However, 
as  the  mechanisms  were  not  damped,  considerable  fluctuations  and  shifts  in  the  ratings  were  noted 
owing  to  pulsating  flow  into  and  out  of  the  float  chambers.  In  particular  the  development  of  standing 
waves  of  varying  intensity  at  differing  stages  around  the  stilling  tubes  meant  that  only  local  water  level 
and  spurious  water  slope  variations  were  indicated.  Instead  water  depth  was  measured  directly  at  half 
meter  intervals  every  half  hour  across  the  section  beneath  the  measuring  platform.  Meaningful  depth- 
discharge  relationships  were  then  devised  for  individual  hydrographs  by  integrating  depth-averaged 
velocity  across  the  section  in  relation  to  the  actual  wetted  area. 

ii)  Velocity:  The  downstream  component  of  velocity  in  the  vertical  was  measured  primarily  using  an 
array  of  six  OTT  c2  velocity  meters.  Vertical  spacings  could  be  altered  with  a  minimum  of  50mm  and 
the  voltage  output  from  each  instrument  was  logged  simultaneously  and  converted  to  calibrated  velocity 
on  a  Hunter  micro-computer.  Velocity  profiles  were  taken  at  0.5m  intervals  across  the  channel  from  the 
measuring  platform  every  hour. 

A  single  head  OTT  electromagnetic  miniature  current  meter  was  also  used  to  obtain  detailed  velocity 
readings  at  centimetre  intervals  in  the  vertical.  This  instrument  has  the  advantage  of  determining  the 
velocity  closer  to  the  bed  than  the  c2  impeller  and  senses  the  flow  over  a  surface  rather  than  integrating 
over  the  area  of  a  rotor.  This  meter  could  be  moved  in  the  vertical  without  altering  the  positioning  of 
the  wading  rod  with  a  precision  of  1mm.  Data  were  stored  directly  on  a  Hunter  micro-computer  or 
latterly  on  a  Toshiba  T3100SX  PC  for  immediate  pre-processing  in  the  field. 

The  turbulence  measuring  system  employed  was  supplied  by  VALEPORT  Instruments,  Dartmouth,  UK, 
and  consisted  of  a  pair  of  17cm  annular  open  bead  EMCM’s  and  associated  electronics.  Using  a  "Y" 
shape  mounting  spar,  the  coil  axis  of  each  EMCM  head  is  fixed  at  45°  to  the  principal  flow  streamline 
and  at  90°  to  the  other  (Fig.  21a).  Prior  to  fieldwork,  calibration  was  undertaken  in  a  large  recirculating 
flume  in  the  Mechanical  Engineering  Department  at  Liverpool  University,  UK.  In  all  cases,  EMCM 
sensitivity  was  approximately  the  same  (ie.,  1  Volt  m'1  s'1),  the  relationship  between  EMCM  output  and 
flow  velocity  was  linear  and  the  standard  deviation  of  the  measurements  at  a  given  velocity  was  very 
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low  (Fig.  21b). 

Measurements  of  orthogonal  turbulent  flow  components  (u,  v  and  w)  were  made  using  the  EMCM 
positioned  at  24cm  above  the  bed.  A  sensitive  Hydrophone  was  attached  close  to  the  current  meter 
heads  on  the  EMCM  spar  at  approximately  23cm  above  the  bed  (Fig.  21a).  Following  laboratory  tests, 
unwanted  EMCM  "noise"  was  reduced  by  passing  output  signals  through  a  15Hz  low-pass  filter.  A 
smoothing  envelope  placed  around  the  output  from  the  log  amplifier  effectively  increased  the  signal  to 
noise  ratio  of  the  hydrophone  without  compromising  the  sensitivity  to  SGN. 

Data  from  the  Gve  output  channels  were  sampled  at  30Hz  using  a  data  logging  system  and  software 
developed  at  the  Proudman  Oceanographic  Laboratory,  UK,  and  stored  directly  to  hard  disk  on  a 
portable  PC.  Following  each  measurement  campaign,  raw  data  were  screened  for  unwanted  noise  spikes. 
In  all  cases,  data  quality  was  high  and  little  editing  was  required.  Figure  22  illustrates  schematically  the 
instrumentation  and  data  acquisition  system  used. 

Measurements  of  flow  turbulence  and  SGN  were  conducted  during  the  period  4  to  8  June  1992.  As 
stage  remained  low  during  this  period  it  was  necessary  to  deepen  artificially  the  stream  channel  and  to 
divert  all  available  flow  through  the  modified  section.  Following  an  adjustment  period  of  several  hours, 
the  bed  appeared  to  re-armour,  relatively  natural  conditions  were  re-established  and  measurements  could 
proceed  in  a  flow  depth  of  approximately  150cm. 


Using  a  purpose-built  mounting  bracket  securely  attached  to  the  large  timbers  of  the  bridge  spanning  the 
upstream  portion  of  the  test  site  (Fig.  5)  it  was  possible  to  position  the  EMCM  and  hydrophone  system 
at  a  fixed  and  stable  position  above  the  irregular  stream  bed  (Fig.  21a).  Despite  efforts  to  ensure  the 
system  was  rigid,  slight  vibration  was  noted  on  the  EMCM  mounting  spar  due  to  sensor  and  support 
oscillation  in  the  fast  flow.  Possible  leakage  of  this  motion  into  EMCM  signals  is  discussed  below.  As 
this  campaign  is  the  first  to  deploy  this  assemblage  of  instrumentation  the  full  details  of  signal 
preprocessing  are  given  in  Annex  1. 

iii)  Measurement  of  bed  elevation  and  water  surface  elevation:  Variation  in  the  bed  and  water  surface 
elevations  were  recorded  using  a  system  of  vertical  rods  spaced  at  0.10m  intervals  across  the  measuring 
platform.  This  apparatus  allows  a  manual  traverse  of  the  stream  taking  readings  of  the  bed  and  water 
surface  elevation  relative  to  a  fixed  datum  at  half  or  hourly  intervals.  The  sampling  error  was 
determined  to  be  +/-  1.0cm  (  De  Jong  pets  comm,  1992). 

d)  Processing  of  hydraulic  data. 

For  brevity  here,  the  preprocessing  of  the  turbulence  records  is  detailed  in  Annex  1.  However,  the  near¬ 
bed  shear  stress  (t),  was  calculated  using  the  Reynolds  stress  (RS)  approximation: 


Page  17  of  46  Pages 


T  ■  -pKW 

or  obtained  using  the  Turbulent  Kinetic  Energy  (TKE)  approximation: 

t  -  0.19p E 


DAJ A45 -90-C-0006 
FINAL  REPORT 


(1) 
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where  p  is  the  density  of  water  and  E  is  defined  in  Annex  1. 
The  bed  shear  velocity  (u.)  is  given  by: 


*.  *  (t/P) 


,0J 


(3) 


Roughness  length  (z„)  and  drag  coefficient  (Q)  were  calculated  from  the  following  approximations: 


<*> 


and; 


where  k  is  von  Kirmin’s  constant  (0.40). 


For  those  near-bed  profiles  which  were  logarithmic  (see  below)  the  near-bed  shear  velocity  and 
roughness  length  were  also  calculated  from  the  equation: 
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where  u,  is  the  velocity  at  a  height  z  above  the  bed,  z,  is  the  reference  height  (typically  0.7DM  )  and  the 
toughness  length  is  the  intercept  of  the  regression  function. 

The  approximate  relationship  between  the  roughness  length  obtained  from  equation  6,  Nikuradse’s 
equivalent  roughness  and  the  grain-size  of  the  bed  material  was  considered  using  the  general  scaling 
relationship; 


3 Qzr  k,  - 


O) 


In  many  gravel-bed  rivers  the  constant  a  is  often  considered  to  be  equal  to  3.551  when  the  representative 
grain-size  (Ds)  is  set  equal  to  DM. 

The  local  shear  stress  (t  =  pu.2)  was  determined  from  individual  profile  data  and  the  Darcy- Weisbach 
hydraulic  roughness  coefficient  (j)  from: 

'-■tel 


where  U  is  the  velocity  averaged  over  the  boundary-layer.  Alternatively,  section-averaged  values  of  all 
parameters  (eg.  x')  were  obtained  from  data  integrated  across  the  section  as  appropriate. 

Total  streampower  (Q)  and  the  energy  slope  (5)  of  the  study  teach  were  calculated  from: 


Q  ■  Ux1  ■  y  QS 


where  y  &  the  specific  weight  of  water  and  Q  is  the  discharge. 
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RESULTS 


Unsteady  Transport 


Spectral  Analysis 

Preliminary  consideration  of  the  original  detector-log  signals  showed  that  at  time-scales  <  ca  300 
seconds  no  temporal  structure  in  the  disposition  of  signals  was  evident  It  follows  that  there  can  be  no 
correlation  to  bulk  flow  parameters  at  this  scale.  The  lack  of  correlation  to  bulk  flow,  and  scale 
considerations,  imply  that  individual  particle  movement  at  the  sub-minute  scale  is  related  to  detailed 
turbulent  flow  structure90,  as  well  as  to  the  vagaries  of  individual  particle  entrainment  and  the  random 
break-up  of  discrete  dusters1142.  At  scales  greater  than  a  few  minutes  (Fig.  7),  apparent  temporally- 
co  be  rent  structures  can  be  detected  visually  although  stochastic  inputs  are  still  strongly  in  evidence. 
Manipulation  of  such  data  is  appropriately  viewed  as  a  problem  in  time  series  analysis.  Sediment 
transport  data  cannot  be  expected  to  show  strong  periodicity5*4443,  in  part  because  the  duration  of  the 
hydrograph  limits  record  length.  Nevertheless  the  raw  spectra,  and  a  series  of  smoothed  variance  spectra, 
were  calculated  for  the  records  from  each  detector  block,  using  the  Parzcn  window  and  10%  taper36. 

Also  calculated  were  the  cross-correlations  between  series  although  not  all  details  are  reported  here.  A 
typical  spectrum  is  shown  in  Fig.  8.  A  number  of  observations  may  be  made.  Firstly,  all  locations  across 
the  channel  width  bad  broadly  similar  spectra  indicating  a  similar  dynamic  response  irrespective  of 
position  across  the  channel;  this  was  confirmed  by  high  cross-correlations  and  near  negligible  phase  shift 
between  series.  Secondly  the  spectra  are  somewhat  flat  indicating  that  variance  is  spread  relatively 
evenly  across  all  frequencies.  An  occasional  significant  peak  may  occur  at  an  approximately  half-hour 
period  but  the  confidence  level  is  not  high.  No  statistically  significant  peaks  occur  in  the  records  at  the 
high  frequency  end  of  the  spectra.  A  significant  peak  at  the  low  frequency  end  of  all  spectra  (<0.01)  is 
associated  with  the  diurnal  periodicity  of  the  hydrographs  as  has  been  noted  before30.  Further  spectral- 
analysis  is  not  warranted,  being  impeded  by  die  length  of  record.  However,  it  is  possible  that  particle- 
count  records  with  a  longer  time-base,  might  demonstrate  other  significant  frequencies.  Bunte®  for 
example  detected  a  frequency  at  0.033  (5  hour  period)  from  visual  inspection  of  a  bedload  record  from 
Squaw  Creek  but  this  was  not  tested  statistically.  Nonetheless,  in  the  present  case  we  can  conclude  that 
no  periodic  pulsing  occurred  during  the  sampled  events;  rather  after  excluding  die  significant  correlation 
with  the  period  of  the  hydrograph,  unsteadiness  in  transport  was  statistically  random.  However,  the 
similarity  of  spectra,  the  strong  cross -correlation  between  locations  and  the  negligible  phase  shift 
indicates  that  the  dynamic  response  of  the  bar-top  environment  essentially  is  contemporaneous  with  the 
talweg  locations.  Notably,  bar-top  dynamics  are  not  driven  by  sediment  moving  slowly  from  die  channel 
onto  the  bar  top.  This  is  an  important  observation  in-as-much  as  it  indicates  a  forcing  function  operative 
at  span-widths  equivalent  to  the  channel  breadth  and  reach  length. 


I 
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Statistical  Distribution  Of  Transport  Rate 

The  final  proposition  in  the  last  section  was  further  tested  by  considering  the  statistical  distribution  of 
the  ten  minute  sampling-period  data.  A  number  of  laboratory  data  sets,  representing  bedload  sampled 
over  0.5  to  1  minute  periods,  have  been  shown  to  broadly  compare  with  a  normal  distribution21-37; 


At)  *  -?=•** 


(10) 


The  normalized  transport  rates  (t)  ate  given  as; 


t  • 


o 


(ID 


where  4,  and  qh  ate  the  individual  bed  load  count  rates  and  the  avenge  rate,  and  a  is  the  standard 
deviation. 

In  contrast,  data  sets  for  instantaneous  sampling,  such  as  nw  particle-count  data  should  include  an 
infinitely  large  number  of  zero  or  near-zero  counts.  In  the  latter  case  a  Hamamori3*  distribution  might  be 
appropriate59-®.  Time-avenging  is  usually  an  essential  element  of  sampling  bedload  and  consequently, 
positively-skewed  normal  distributions  should  be  common  with  a  tail  of  large  transport  ntes.  Laboratory 
data  presented  by  Kuhnle  and  Southard21  in  fact  show  weak  positive  skewness  probably  induced  by  the 
sampling  strategy.  The  present  field  data  are  also  well  described  by  skewed  normal  distributions  (Fig. 

9).  The  distributions  when  comparing  both  upstream  and  downstream  records  and  between  events 
broadly  are  similar  in  peakedness  and  degree  of  skew;  implying  comparable  generating  mechanisms  are 
involved.  The  minor  differences  in  the  distributions  ate  more  likely  associated  with  inadequacies  in 
sample  size  than  process  differences.  An  example  is  the  weak  development  of  polymodality  towards  the 
right-hand  side  of  the  distributions.  However  strong  bi-  or  poly-modality  would  be  significant  for  the 
following  reason.  If  distinct  pulsing  occurred  in  the  temporal  distribution  of  time-avenged  bedload 
records,  then  bi  modal  histograms  might  occur.  This  is  because  very  low-rates  and  high  transport  pulses 
would  be  well  represented  whilst  few  data  would  exist  for  medium  transport  ntes.  The  data  from 
individual  sections  across  the  river  are  also  skewed,  but  interestingly,  in  the  case  of  the  23/24th  May 
data,  the  bar  and  interface  data  have  a  greater  weighting  to  smaller  tn  ns  port  events  than  is  observed  for 
the  channel  sections;  this  is  especially  marked  in  the  case  of  the  interface  (Figs.  10,11  &  12).  The  data 
for  the  5/6th  June  demonstrate  changes  in  the  gnin  count  distributions  induced  by  tn  ns  port  between 
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upstream  and  downstream  detectors.  Distributions  become  increasingly  peaked  and  less  skewed; 
especially  over  the  bar-top  (eg  Fig.  12).  Although  substantial  variation  exists  in  the  transport  rates 
across  the  sections,  there  is  no  evidence  that  the  coarse  bedload  moves  in  distinct  bands  with  near  zero- 
transport  in  between  (cf.  61“).  Although  maximum  transport  is  associated  with  the  talweg  location  the 
whole  stream-width  is  active  once  bedload  motion  begins  (Figs.  7  &  13). 

General  Transfer  Function 


Viewing  the  transport  event  as  a  whole,  continuity  in  the  mass  of  sediment  transported  as  a  sheet  may 
be  maintained  such  that  as  many  particles  may  exit  the  reach  as  enter.  The  dynamic  of  the  particles 
within  a  sheet  may  be  seen  as  analogous  to  the  passage  of  a  tracked-vehicle,  whereby  individual  plates 
are  laid-down  at  the  front  subsequently  to  be  picked-up  again  at  the  rear.  A  comparison  of  the  upstream 
particle  count  (17479)  for  the  whole  event  with  the  downstream  count  (16319)  can  be  used  as  a  measure 
of  continuity  and  indicates  that  some  deposition  occurs  within  the  reach.  However,  in  this  example,  the 
difference  of  only  7%  could  be  owing  to  instrument  error. 

If  continuity  of  transport  is  maintained  over  the  short  study  reach,  it  follows  that  the  upstream  count 
series  should  be  capable  of  generating  the  downstream  series.  Visual  inspection  of  the  stream-wide 
series  for  the  23/24th  event  shows  an  initial  in-phase  correlation  of  transport  peaks  which  later  become 
lagged  with  the  downstream  detector  peaks  occurring  up  to  ten  minutes  later.  By  treating  both  upstream 
and  downstream  data  series  as  an  input-output  model  using  recursive  iteration63,64,  in  principle,  general 
sediment  transfer  functions  can  be  defined  wherein  stable  model  structure  might  elucidate  the  physical 
nature  of  sheet  dynamics.  The  statistical  analysis  was  implemented  using  a  computer  package  called 
microCAPTAIN65,  in  an  attempt  to  find  a  black-box  transfer  function  capable  of  generating  the 
downstream  record  from  the  upstream  series.  The  results  for  the  whole  event  were  poor,  mainly  because 
for  this  particular  date  set  it  proved  difficult  to  mode]  the  initial  spike  in  both  records  (Fig.  14). 
However,  consideration  of  sub-sets  (Fig.  IS)  proved  more  successful  (Fig.  16).  For  example,  modelling 
from  the  SSth  to  99  points  (ie  from  1638hr  to  2302hr),  where  a  lag  of  approximately  one  time-step  is 
evident,  gave  a  significant  coefficient  of  determination  of  0.61  (Young  Information  Content  a  -4.627).  A 
useful  introduction  to  the  statistical  treatment  of  similar  date  sets  is  given  by  Young  and  Wallis66.  The 
appropriate  function  in  this  example  is: 

/*-!.!  1583/t_,40.4750i*_2*0.2504«*_1  (12) 


where 
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0.2504ut_, 

1 -1.1 1583Z'1  *0.47503z -2 


(13) 


Where  ut.u  is  the  lagged  input,  lk  is  the  prediction  at  time  k  and  z  is  an  operator. 

With  only  the  one  data  set  available  at  this  time  it  is  not  possible  to  say  whether  the  behaviour  noted  is 
characteristic.  Clearly  the  lag  identified  is  dictated  in  part  by  the  averaging  period  used  in  pre-processing 
the  count  data.  However,  the  procedures  adopted  above  serve  to  establish  a  framework  within  which 
similar  records  usefully  can  be  examined  once  these  become  available.  The  coherent  structure  of  the 
dynamic  sheet  moving  through  the  reach  does  tend  to  indicate,  however,  that  there  is  little  significant 
exchange  of  material  between  the  sheet  and  the  bed  material.  This  observation  is  explored  more  fully 
below. 

Bed  Level  Dynamics 

The  movement  of  bedload  sheets  through  the  study  reach  may  be  reflected  in  bed  level  fluctuations  if 
exchange  of  sediment  occurs  in  the  vertical.  Otherwise  the  relationship  of  the  sheet  and  the  bed  would 
be  passive,  with  bedload  transport  taking  place  over  a  static  surface.  Analysis  of  the  topographic  data 
showed  that  slight  fluctuations  in  bed  level  did  occur  so  that  the  sheets  interacted  with  the  bed  material. 
Such  a  situation  is  superficially  akin  to  the  Phase  /  transport  described  by  Jackson  and  Bescbta3,  where 
fine  sediment  is  entrained  from  and  over-passes  an  intact  armour  layer  and  only  local  readjustment  takes 
place,  but  differs  in-as-much  as  here  coarse  bedload  (up  to  127mm)  over-passes  a  largely  static  armour 
and,  as  will  be  shown  below,  fine  sediment  is  deposited  into  the  armour. 

At  each  sample  point,  the  minimum  bed-level  elevations  for  each  event  were  subtracted  from  the 
maxima  (Fig.  17a)  to  provide  a  measure  of  the  thickness  of  the  active  layer.  Interestingly  the  coarse- 
bedded  talweg  changed  little  whilst  the  bar  top  changed  only  markedly  in  the  vicinity  of  the  2m  mark 
where  a  slight  topographic  low  already  existed  at  the  beginning  of  the  event  Maximum  change  in 
elevation  occurred  along  the  interface  between  the  bar  ami  the  talweg  channel  (Fig.  17b)  and  is 
equivalent  to  the  disturbance  of  one  or  two  armour-particle  diameters.  Elsewhere  changes  are 
attributable  to  isolated  exposed  particles  being  mobilized.  The  spatial  distribution  of  these  changes 
through  time  are  illustrated  in  Fig.  18.  However,  following  Hoey  and  Sutherland2,  the  main  temporal 
changes  are  illustrated  well  by  considering  the  standard  deviation  of  the  bed-level  records  (Figs.  19  & 
20).  In  both  examples,  the  interface  area  is  more  dynamic  than  bar-top  or  channel  talweg,  but  it  is 
notable  that  as  the  interface  roughens,  the  bar  top  is  smoothed  and  vice  versa.  The  variance  in  the  data 
declines  slightly  as  bedload  transport  commences  and  is  then  especially  great  in  the  period  following  the 
passage  of  a  bedload  sheet  through  the  reach  (Fig.  20).  Consequently,  the  bed  becomes  rougher  once 
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bedload  transport  ceases.  Such  physical  adjustments  are  explained  by  the  initial  infilling  of  the 
interstices  of  a  largely  static  cobble  surface  with  finer  pebble-sized  material;  material  in  particular  being 
exchanged  between  interface  and  the  bar-top  area.  As  noted  above,  this  fine  bedload  is  derived  from  up¬ 
stream  of  the  study  reach  and  transported  into  the  reach.  Subsequently,  once  the  peak  of  the  bedload 
wave  passes,  the  flow  is  still  competent  to  transport  more  particles  than  it  now  contains  and  the  recently 
deposited  material  is  reentrained  resulting  in  a  coarsening  of  the  bed.  The  infilling  process  however  is 
very  dynamic  in-as-much-as  many  individual  pebbles  are  probably  deposited,  reentrained,  and  then 
replaced  by  others  constantly.  The  temporal  changes  in  bed  roughness  will  have  an  effect  on  turbulence 
intensity  over  the  bar-top  and  concomitant  changes  should  be  expected  in  bulk  flow  parameters.  The 
latter  proposition  is  explored  below. 

Hydraulics 

Turbulence 


Statistical  Aspects: 

In  Figures  23  to  27,  for  clarity,  time  series  plots  for  the  first  0.55  minutes  of  runs  4  to  8  show  speed,  u, 
v,  w,  uw,  vw,  stress  magnitude  (x),  stress  direction  (W)  and  SGN.  In  addition  power  spectre  and 
cospectra  obtained  using  8192  data  values  for  u,  v,  w,  uw  and  vw  are  also  shown.  A  summary  of 
average  current  speed,  RMS  turbulence  intensity  for  u,  v,  and  w,  bed  shear  velocity  obtained  using 
turbulent  kinetic  energy  (TKE)  and  Reynolds  stress  (RS)  approximations,  dreg  coefficient  (Cd)  and 
roughness  length  (zn)  values  are  given  in  Table  2. 

» 

As  expected,  the  statistical  distribution  of  instantaneous  u,  v  and  w  values  are  found  to  be  Gaussian 
(Fig.  28).  In  contrast,  uw  and  vw  and  SGN  series  are  found  to  be  positively  skewed  and  are 
characterized  by  a  long  tailed  distribution  (Fig.  29a  &  b).  Although  not  demonstrated  here,  previous 
studies  indicate  that  the  high  stress  values  associated  with  infrequent  "events”  shown  in  Fig.  29a  are 
likely  to  be  positively  correlated  with  sediment  mobilization34.  Further  discussion  and  statistical  analysis 
of  these  distributions  is  given  by  Williams  and  Tawn30.  The  high  frequency,  low  magnitude  SGN  events 
illustrated  in  Fig.  29b  are  considered  to  be  associated  primarily  with  flow  'noise*  rather  than  with 
sediment  transport.  In  contrast,  high  magnitude,  low  frequency  events  are  considered  to  originate 
through  grain  impacts  or  other  acoustic  sources  at  distances  away  from  the  hydrophone  and  show  little 
correlation  with  local  hydrodynamic  conditions.  Further  consideration  of  these  observations  is  given 
below. 

The  standard  deviation  for  individual  10  minute  blocks  of  EMCM  data  from  a  continuous  measurement 
period  of  1.5  hours  showed  no  significant  variation  (<  1cm  s'1)  through  time.  Further,  although  spectra 
obtained  for  the  whole  record  indicated  small  rises  or  declines  in  stream  stage,  they  did  not  exhibit  low 
frequency  energy  peaks  (5-10  minutes)  as  noted  in  other  geophysical  flows676*.  In  general,  as  the  typical 
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turbulent  eddy  sizes  scale  with  flow  depth  (see  below),  fluid  motions  with  a  period  of  this  order 
probably  originate  through  a  mechanism  unassociated  with  normal  production  and  dissipation  processes 
for  turbulent  eddies.  Whilst  not  yet  confirmed,  it  has  been  suggested  by  Lapointe®8  and  others  that  low 
frequency  turbulent  motions  in  marine,  estuarine  and  large  river  situations  may  be  associated  in  some 
way  with  bedform  interactions  with  flow  during  active  sediment  transport.  Although  differing  in  flow 
Reynolds  number  by  a  factor  of  approximately  10  from  the  shallow  marine  situation  (z  =  20m;  U  =  1  m 
s'1)  and  with  turbulent  eddies  scaling  more  appropriately  with  flow  depth  than  with  grain  size  (see 
below),  failure  to  detect  these  motions  in  Squaw  Creek,  where  sediment  transportation  was  negligible 
and  recognizable  bedforms  were  absent,  would  appear  to  support  this  suggestion. 

In  common  with  u,  v  and  w  time  series,  instantaneous  velocity  changes  recorded  by  the  EMCM  system 
are  found  to  be  Gaussian  and  indicate  that  particles  on  the  bed  are  likely  to  be  subjected  to  rapid,  high 
magnitude  flow  accelerations  and  de-accelerations  through  time  (Fig.  30).  DeGning  the  period  of  an 
"event”  as  being  the  time  during  which  the  velocity  of  a  given  flow  component  continues  to  rise  or  fall, 
Fig.  31  shows  that  short  duration  events  (<2.5  sec)  give  rise  to  the  largest  acceleration/de-acceleration 
terms  and  exhibit  an  approximately  exponential  decay  with  event  length. 

When  normalizing  using  U.^  values  (Table  2)  RMS  turbulence  intensity  values  show  a  high  degree  of 
consistency  and  range  between  1.9  -  2.3,  2.2  -  23  and  1.0  -  1.6  for  u,  v  and  w  respectively.  These 
values  are  in  good  agreement  with  those  obtained  in  other  geophysical  flows69.  RMS  turbulence  intensity 
values  normalized  using  U>7XE,  however,  are  consistently  lower.  Relatively  steady  flows  during  the 
measurement  periods  prevented  study  of  the  RMS  turbulence  levels  over  a  range  of  flow  contritions. 

As  time  series  obtained  in  this  study  correspond  more  closely  to  spatial  than  temporal  variations  in 
velocity14,  power  spectra  obtained  using  a  fast  Fourier  transform  (FFT)  were  converted  to  wavenumber 
spectra  (not  illustrated)  by  the  transform: 


E(.k)  «  (U(z)/2n}S(fl 


(14) 


where  E(k)  is  the  wavenumber  spectrum,  k  =  2jtf/U(z),  z  is  the  sensor  height  above  the  bed  and  S(f)  is 
the  power  spectrum  at  frequency  f.  As  frequency  is  the  reciprocal  of  period  (T),  the  resulting  spectra 
plotted  in  a  normalized  form  can  be  used  to  estimate  typical  eddy  length  scale  L  by; 


*  -  In/UT  -  2n/t 


(15) 
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Thus 


L  =  Inxjk 


(16) 


At  peak  spectral  values,  typical  turbulent  eddy  length  scales  for  u,  v  and  w  were  found  to  be 
approximately  0.52,  0.54,  and  0.11m  respectively.  Similar  results  were  also  obtained  through 
autocorrelation  analysis  of  the  time  series,  (not  illustrated  here),  and  indicate  clearly  that  typical  eddies 
scale  with  flow  depth  in  Squaw  Creek. 

Spectra  plotted  in  Figures  23  to  27  exhibit  spectral  slope  values  ranging  between  -1.85  and  2.58  and 
deviate  from  the  expected  value  of  -5/3.  These  results  demonstrate  that  fully  isotropic  conditions  do  not 
occur  at  the  maximum  spectral  wavenumbers  in  this  study.  These  discrepancies  are  considered  to  arise 
owing  to  the  physical  size  of  the  EMCM  heads  and  to  sampling  interval  whereby  the  small  scale 
turbulent  motions  responsible  for  energy  dissipation  are  not  recorded.  This  is  a  common  problem  with 
all  current  measurement  devices  employed  in  field  studies  of  turbulence  and  can  lead  to  underestimation 
of  shear  stresses  when  employing  TKE  and  RS  methods70. 

Evident  on  all  spectra  obtained  in  Squaw  Creek  are  small  peaks  at  frequency  values  of  approximately  5 
Hz.  As  these  are  absent  in  calibration  spectra  their  presence  may  be  attributable  to  the  combined  offset 
of  EMCM  head  vibration  in  the  relatively  fast  flow  and  to  regular  eddy  shedding  phenomena  upstream 
of  the  test  section.  Although  subject  to  rapid  oscillation  in  the  calibration  flume,  the  narrow  frequency 
band  associated  with  the  present  spikes,  however,  suggests  vibration  of  EMCM’s  is  the  primary  cause. 
At  worst,  presence  of  these  spurious  "turbulence”  signals  results  in  an  over  estimation  of  bed  shear 
stress  of  <  2%  and  can  be  ignored  for  the  purposes  of  the  present  study. 

Hydraulic  parameters  derived  from  turbulence  data: 

It  is  generally  accepted  that  the  "constant  stress  layer"  of  a  turbulent  boundary  layer  over  an 
hydraulically  smooth  surface  occupies  approximately  the  first  10-15%  of  the  total  layer  thickness  and 
measurements  of  Reynolds  stresses  within  this  region  yield  an  approximation  to  bed  shear  stress.  Given 
a  flow  depth  at  the  present  study  site  of  approximately  150cm,  Reynolds  stress  approximations  are  only 
likely  to  be  valid  for  z  up  to  15cm  if  the  conditions  of  hydraulic  smoothness  can  be  met.  Clearly  given 
the  nature  and  arrangement  of  the  coarse  particles  comprising  the  bed,  this  assumption  is  likely  to  be 
misleading  owing  to  the  presence  of  an  internal  boundary  layer  with  turbulence  characteristics  probably 
scaling  with  a  bed  toughness  parameter.  It  would  not  be  expected,  therefore,  that  the  present 
measurements  obtained  at  z  *  24cm  using  relatively  large  EMCM  beads  will  give  results  that  relate  to 
bed  shear  stresses  directly.  In  combination  with  sampling  errors  associated  with  the  present  EMCM 
system,  it  would  not  be  surprising,  therefore  to  find  that  u.  values  in  Table  2  differ  from  those  obtained 
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using  a  log-profile  fitting  technique. 

Values  of  u.  in  Table  2  obtained  using  TKE  and  RS  methods,  however,  are  not  only  very  consistent 
during  measurement  runs  4  to  8,  but  also  agree  well  with  some  u.  values  obtained  using  the  log-profile 
(Fig.  36).  On  the  basis  of  these  results  it  would  appear  that  the  present  EMCM  system  provides  a 
realistic  estimation  of  bed  shear  stress  in  highly  turbulent  flow  conditions  encountered  in  Squaw  Creek. 

Also  shown  in  Table  2  are  drag  coefficient  (C^)  and  roughness  length  (zj  values  obtained  using 
equations  4  and  S.  In  common  with  u.  values,  Cd  and  values  are  generally  consistent  over  runs  4  to  8 
and  are  generally  lower  than  would  normally  be  expected  given  the  grain  size  distribution  of  the  bed 
sediments.  In  contrast  zD  values  obtained  from  velocity  profile  data  (reported  below)  are  consistent  with 
expectations  based  on  a  consideration  of  bed  grain  size  and  equation  7.  The  unusually  low  values  of  z0 
obtained  from  the  turbulence  data  series  will  require  further  investigation.  One  possible  explanation  is 
that  in  high  velocity  flow  over  rough  beds,  skimming  flow  develops  such  that  roughness  lengths 
obtained  using  outer-layer  turbulence  measurements  are  related  to  the  interface  between  the  outer-layer 
and  the  near-bed  boundary-layer.  Consequently,  no  direct  relationship  should  then  be  expected  between 
z„  and  the  bed  material  characteristics. 

SGN  measurements: 

Time  series  plots  of  SGN  in  Figs.  23  to  27  show  little  visual  correspondence  with  any  flow  parameter 
obtained  using  the  EMCM  measurements.  This  has  been  confirmed  from  cross-correlation  analyses  using 
30  minute  time  series.  Further,  the  spectrum  of  SGN  is  closely  similar  to  that  for  u  (Fig.  32),  and  for 
SGN  output  in  flows  without  sediment  transport.  This  suggests  strongly  that  the  hydrophone  is  only 
recording  flow  "noise".  Attempts  to  seed  Squaw  Creek  with  copious  amounts  of  mixed  grain  size 
particles  on  a  number  of  experimental  runs  failed  to  produce  correlation  between  SGN  and  bedload 
transport  (Fig.  32).  These  results  were  recognized  in  the  field  during  preliminary  data  analysis  and  thus 
an  alternative  experimental  technique  was  sought. 

Selected  particles  from  the  bed  were  collected  and  sorted  by  size.  By  gently  knocking  particles  together 
underwater  at  a  distance  of  0.5m  away  from  the  hydrophone  it  was  possible  to  detect  and  record  SGN 
signal  for  a  range  of  sizes  and  positions  from  the  sensor.  During  experimental  runs  over  approximately 
one  minute,  attempts  were  made  to  reproduce  regular  collisions  between  the  test  particles  in  ordn  that 
they  may  be  more  clearly  identified  in  the  hydrophone  record.  In  the  event  this  proved  to  be 
unnecessary  as  the  resulting  SGN  level  was  found  to  be  significantly  above  the  background  acoustic 
noise  of  Squaw  Creek. 

SGN  data  typifying  an  impact  test  at  1.0m  from  the  hydrophone  using  well-  rounded  particles  with  a 
diameter  in  the  range  32  to  64mm  is  shown  in  Fig.  33.  In  common  with  previous  plots,  SGN  is  shown 
on  a  logarithmic  scale.  The  near  regular  peaks  present  in  Fig.  33  rise  by  at  least  two  orders  of 
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magnitude  above  the  "background"  level  and  thereby  provide  convincing  evidence  that  SGN  can  be 
detected  in  Squaw  Creek.  This  is  illustrated  further  in  Fig.  34  which  shows  SGN  and  u  spectra  on  the 
same  axes.  Here,  in  contrast  with  Fig.  32,  the  spectral  characteristics  of  u  and  SGN  differ  significantly 
over  the  range  of  frequencies  shown.  Given  these  clear  results  and  the  lack  of  natural  bedload  transport, 
no  further  analysis  of  SGN  data  was  attempted. 

Depth-integrated  parameters 

Shallow  flow  over  large  bed  roughness  may  be  associated  with  local  flow  accelerations  and  velocity 
profile  distortion  such  that  uniform  flow  assumptions  are  violated.  The  effect  of  large  relative-depth 
ratios  on  the  velocity  distribution  throughout  the  section  is  examined  for  data  sets  obtained  at  half- 
hourly  intervals  throughout  the  hydrographs  (eg.  Fig.  35).  Water  surface  elevation  closely  matches 
variation  in  the  topography  across  the  section  and  isovels  are  distorted  and  close  together  in  the  vicinity 
of  topographic  high  points.  These  high  points  were  often  isolated  blocks  around  and  over  which  the  flow 
is  forced  to  deviate.  Although  this  induces  strong  secondary  currents  in  the  immediate  vicinity  to  the 
blocks,  there  was  no  indication  that  such  coherent  secondary  flow  cells  were  persistent  stream-wise. 

Such  effects  tended  to  reduce  as  water  depth  increased  and  the  flow  became  more  two-dimensional  with 
a  water  surface  less  distorted.  Variations  on  this  theme  occurred  owing  to  the  local  development  of 
standing  waves  out-of-phase  with  the  bed  topography,  but  owing  to  low  Froude  numbers,  this  was 
exceptional. 

It  has  been  established  that  if  strongly  local  flow  effects  are  avoided,  then  with  care,  results  obtained 
from  log-profile  analysis  can  be  representative  of  the  bulk  flow  hydraulics71,72.  Before  calculating 
parameter  values  from  velocity  profile  data  obtained  using  the  array  of  six  impellers,  a  careful 
examination  of  the  vertical  flow  structure  was  undertaken  using  a  vertical  deployment  of  the 
electromagnetic  sensor.  With  this  instrument  velocity  was  measured  at  half  to  one  centimetre  intervals  in 
the  vertical  (Fig.  36)  with  the  lowest  reading  taken  by  nestling  the  sensor  between  surface  bed  particles 
such  that  the  zero  datum  was  ca  0.7DW.  A  number  of  observations  may  be  made.  The  profile  extending 
to  ca  25%  of  the  depth  often  is  close  to  logarithmic.  The  point  of  inflection  from  this  log-profile  is 
representative  of  the  thickness  of  the  constant  stress  boundary-layer  associated  with  the  grain  roughness, 
whilst  the  velocity  readings  above  this  point  are  representative  of  an  outer-layer  associated  with  both 
grain  and  form  roughness.  In  this  outer  layer  the  vertical  velocity  structure  demonstrated  greater 
variance  about  respective  trend-lines  than  was  noted  within  the  boundary-layer.  Following  Wilkinson73, 
only  those  points  representative  of  the  lower  15%  of  the  depth  and  within  a  few  centimetres  of  the  bed 
were  used  to  calculate  near-bed  parameters.  This  procedure  generally  gave  profiles  close  to  logarithmic 
with  little  data  scatter  about  the  regression-line  (r  values  >0.95).  Following  Jackson74  consideration  was 
given  to  adjusting  the  zero  datum,  but  this  was  usually  found  to  be  inappropriate  because  the  goodness- 
of-fit  of  data  to  a  log-profile  was  excellent  without  any  adjustment  and  it  was  concluded  that  the  zero- 
plane  was  well  represented  by  the  mode  of  deployment  (see  also  Soulsby®).  There  is  no  objective 
method  available  for  making  an  adjustment  in  shallow  flow  over  coarse  gravel  and  instead,  those  near- 
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bed  profiles  (eg  Profile  7)  which  did  not  meet  the  log-normal  criteria  were  not  analyzed  further. 

Inflected  profiles  are  common  in  the  fluvial  environment  and  often  may  be  related  to  flow  over  surface- 
roughness  changes73,  where  concave-up  profiles  reflect  a  rapid  transition  in  the  development  of  the 
thickness  of  the  boundary-layer  as  flow  traverses  from  smooth  to  rough  beds  and  vice  versed.  Of  the 
non-logarithmic  distributions,  convex-up  profiles  were  more  common  than  concave-up  profiles  but  no 
evident  downstream  change  in  surface  roughness  was  present.  In  addition  it  was  noted  that  the  log-layer 
is  often  well  developed  to  circa  25%  of  the  depth  which  mitigates  against  an  explanation  based  on 
invoking  rapidly  reforming  boundary-layers  in  response  to  stream-wise  variation  in  roughness.  Instead, 
these  profiles  are  superficially  similar  to  those  compound  profiles  described  for  developed  flow  over 
coarse  gravel76,77,78  and  by  Dyer79  ,  amongst  others,  for  flow  over  dune  fields.  Consequently,  such 
profiles  can  be  attributed  to  flow  over  composite  roughness;  ie  grain  roughness  and  additional  form 
roughness,  where  the  latter  is  owing  to  isolated  large  blocks80,  or  other  large  bedforms  where  the 
characteristic  outer-region  roughness  lengths  are  greater  than  the  grain  roughness81.  A  particularly 
valuable  mathematical  derivation  of  the  velocity  profile  in  this  environment  is  given  by  Jackson74. 
Immediately  above  the  reference  plane  the  flow  structure  depends  primarily  upon  grain  roughness  and 
only  weakly  upon  the  geometry  of  large-scale  roughness  elements.  In  the  outer  layer,  in  contrast,  the 
flow  structure  depends  to  a  greater  part  on  the  presence  of  large-scale  roughness.  Visual  inspection  of 
the  bed  in  Squaw  Creek  at  the  time  indicated  that  the  compound  profiles  could  be  ascribed  to  wake 
effects  induced  by  larger  than  average  bed-elements  distributed  randomly  across  the  section.  The 
sheltering  effect  of  upstream  blocks  lead  to  a  flow  structure  near  the  gravel-bed  distinct  from  that  above 
the  block  height.  The  values  of  z„  obtained  from  profiles  taken  within  a  few  centimetres  of  the  bed  (Fig. 
36)  are  typically  of  an  order  of  0.55cm;  which  from  Equation  7  indicates  that  k,  for  the  cobble  bed  is 
circa  16cm  or  1.8DW  of  the  bar-top  gravels.  This  value  for  k,  compares  favourably  with  experimental 
data82  obtained  over  similar-sized  cobbles  in  a  laboratory  flume.  Profile  data  taken  at  greater  heights 
above  the  bed,  which  may  include  outer-layer  observations  gave  a  median  zD  value  of  1.6cm  (event  of 
23/24tb  May)  and  1.8cm  (event  of  5/6tb  June)  which  indicates  that  k,  for  the  large  scale  roughness  is  of 
the  order  of  half  a  meter  or  2.5DM  of  the  channel  armour.  It  has  already  been  noted  above  that  these  z„ 
data  obtained  from  log-profiles  are  at  variance  with  the  z0  data  obtained  from  analysis  of  the  turbulence 
record. 

The  vertical  structure  of  air-flow  over  and  through  vegetation  canopies83  also  is  qualitatively  similar  to 
the  form  of  profiles  described  above.  This  analogy  between  flow  over  blocks  on  a  gravel-bed,  dunes  and 
through  vegetation  may  provide  further  insight  into  the  quantitative  description  of  compound  profiles 
within  coarse-gravel  rivers.  It  is  the  authors  intention  to  consider  this  further  when  additional  data  are 
available. 

Interaction  of  bed  material  and  hydraulics 

The  relationship  of  hydraulic  parameters  to  bed  elevation  and  bedload  data  was  considered  both  in  terms 
of  the  spatial  variance  considering  the  individual  data  points  across  the  section  and  for  section-integrated 
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data.  The  variation  in  shear  velocity  for  example  across  the  section  and  through  time  is  shown  in  Fig. 
37.  Correlation  matrices  were  calculated  relating;  (i)  adjacent  bed  elevation  data,  (ii)  adjacent  shear 
stress  data  and,  (iii)  bed  elevation  data  with  shear  stress  data.  These  however  are  not  presented  here  as 
they  showed  little  evidence  of  a  consistent  response  of  the  bed  surface  to  local  variation  in  the  shear 
velocity,  although  there  was  weak  evidence  for  correlation  of  parameters,  and  hence  a  dynamic 
response,  at  the  bar  interface. 

Negligible  entrainment  of  bar  and  channel  gravels  in  the  study  reach  is  further  indicated  by  a 
consideration  of  theoretical  thresholds  for  particle  entrainment.  Carling  and  others82  provide  an 
entrainment  function  for  coarse  ellipsoidal  gravel  in  shallow  streams  such  as  Squaw  Creek: 


2P(P,-P)g  rinft-p)  1 

*Cdn(aJ>2cJ  «*(♦)  [S.751og(3Qz^/Jk^)]2 


For  conciseness  notation  is  given  by  Carling  and  others82  but  here  appropriate  values  are:  Cd  =  1.25,  $  = 
53°,  p  =  2°,  (p.-p)  =  l-65g  cm'3,  g  =  981cm  s'2,  zp=  3.9cm  and  kt  =  15cm.  V  =  4/3jia2 (where  a2  = 
D/Z  the  short  semi-axis  of  the  ellipsoid  and  b  and  c  represent  the  other  major  semi-axes).  Entering  these 
data  into  equation  17  and  solving  for  the  critical  shear  stress  for  entrainment  for  a  variety  of  grain- 
sizes  can  be  ascertained.  The  bed  armour  over  the  bar  top  can  be  characterized  by  a  DM  of  at  least 
90mm  (Table  1).  A  shear  stress  of  circa  75  N  m'2  is  required  to  entrain  any  exposed  armour  particles 
seated  on  shallow  cols  between  coplanar  particles  (exposure  parameter  *P  =  0.84).  Near  coplanar 
particles  (embedded,  VP  =  0.5)  which  form  the  bulk  of  the  armour-layer  require  at  least  150  N  m'2  for 
entrainment.  In  contrast  the  gravel  in  the  sheet  moving  over  the  bar  in  exposed  locations  -  0.99)  can 
be  typified  by  a  DM  of  between  60  and  80mm  (Fig.  38)  which  requires  between  40  and  65  N  m'2  to 
maintain  mobility.  A  sensitivity  analysis,  varying  parameters  between  known  limits,  readily  shows  that 
although  threshold  values  change  slightly,  the  bed  armour  would  remain  intact  (eg  for  the  event  of  23/24 
May  the  measured  shear  stress  varied  between  80  and  130  N  m2  ).  Consequently  all  grain-sizes  within 
the  sheet  should  pass  rapidly  through  the  reach  with  little  entrainment  of  the  compact  bed  armour. 

Section-integrated  data  (Figs.  39  &  40)  showed  conclusively,  that  the  passage  of  bedload  sheets  through 
the  reach  is  indeed  poorly  correlated  to  bed  elevation  but  has  aneffect  on  bulk  hydraulic  parameters. 

This  is  best  illustrated  by  the  event  of  23rd/24th  May  1991.  Although  discharge  tended  to  increase  after 
2000  hrs,  the  arrival  and  passage  of  a  bedload  sheet  through  the  reach  led  to  smoothing  of  the  bed 
surface  in  the  talweg  and  bar-interface  (Figs.  19  &  40)  such  that  velocity  increased,  water  depth  reduced 
and  the  shear  velocity  and  streampower  fell  as  the  hydraulic  roughness  was  reduced  (Fig.  41).  Similar 
reductions  in  bed  roughness  were  associated  with  bedload  transport  during  the  event  of  21st  May  and 
5th/6tfa  of  June  1991  and  reflect  sudden  changes  in  the  energy  slope  as  bedload  sheets  pass  through  the 
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reach.  Following  the  passage  of  the  sheets,  the  bed  physically  and  hydraulically  roughened  (Fig.  20  and 
41)  with  a  concomitant  adjustment  of  other  hydraulic  parameters  including  depth. 

Owing  to  the  coarse  nature  of  the  bedload  and  the  inability  of  the  net  sampler  to  collect  efficiently  the 
fractions  finer  than  11mm,  the  bedload  grain-size  curves  do  not  show  well  the  changes  in  the  finest 
fractions  as  materia]  is  exchanged  with  the  bar  top  (Fig.  38).  The  DM  in  fact  changes  little  (Fig.  42)  but 
there  is  some  indication  that  between  ISOOhr  and  0600hr  the  sheet  at  first  is  dominated  by  finer  gravels 
(1600hr)  with  a  negligible  percentage  coarser  than  84mm.  The  distribution  then  coarsens  slightly  as  it 
enters  the  reach  and  loses  fines  to  the  surface  armour  (1800hr)  after  which  there  is  little  change  in  the 
distribution  until  0200hr,  following  which,  the  tail  of  the  sheet  passes  over  the  bar  and  finer  gravels 
increasingly  dominate  the  distributions.  Consequently,  the  coarse  end  member  declines  steadily  from 
<127mm  to  <45mm  (Figs.  38  &  42)  with  only  sporadic  examples  of  cobbles  in  transport  after  0900hr. 
The  main  sheet  had  passed  by  0800hr  and  transport  rates  dropped  rapidly.  Occasional  bedload  samples 
taken  subsequently  showed  variable  grain-size  distributions.  These  latter  samples  reflect  the  random 
passage  of  a  few  cobbles  and  fines. 

The  variation  in  bed  roughness  and  consequent  changes  in  hydraulic  parameters  during  the  passage  of 
bedload  sheets  demonstrates  that  even  when  the  intensity  of  transport  is  only  low,  and  the  armour  is 
largely  intact,  uniform-flow  resistance  equations  are  unlikely  to  be  applicable.  This  notion  was  tested  by 
applying  Hey’s51  flow  resistance  equation  to  the  data  for  the  23/24th  of  May  and  5/6th  June  when 
transport  was  relatively  weak. 


(1*) 


Hey’s  equation  is  currently  the  most  satisfactory  for  gravel-bed  rivers  giving  acceptable  results  when 
compared  with  measured  flow  resistance  over  stable  beds72.  In  using  equation  18,  it  was  assumed  that 
for  a  broad  shallow  channel51,  a  =  11.1;  whilst  the  fixed  bed  DM  should  be  optimized.  This  latter 
objective  was  achieved  by  first  substituting  observed  values  of  R  and  /  into  the  equation  to  produce  a 
range  of  values  for  DM.  For  example,  the  average  for  the  event  of  5/6tb  June  equalled  197mm  (sd 
33.91mm  for  9  observations).  These  values  are  consistent  with  variability  in  the  armour-layer  grain-size. 
However,  when  (8//)0J  is  calculated  assuming  this  constant  value  of  DM,  for  observed  variation  in  flow 
depth,  the  inadequacy  of  the  function  becomes  apparent  Sediment  transport  commenced  between  the 
first  and  second  observations  smoothing  the  bed  such  that  the  second  observed  value  of  (8/0°^  is  much 
larger  than  predicted  (Fig.  43).  As  bedload  transport  ceased  the  discrepancy  declines,  only  to  begin 
increasing  again  when  transport  starts  again  after  2100hr. 
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4.0  CONCLUDING  DISCUSSION 

■  During  the  relatively  low  flows  encountered  in  1991  and  1992  it  has  proven  possible  to  obtain  good 
quality  measurements  of  flow  turbulence  and  weak  transport  of  bedload  sheets.  Further,  through  the 
local  generation  of  SGN  and  magnetic  induction  it  has  been  shown  that  acoustic  and  magnetic  detection 
of  bedload  in  fast  flowing,  gravel  bed  streams  is  possible  using  the  present  detection  systems. 

■  With  appropriate  scaling,  the  statistical  properties  of  u,  v  and  w  time  series  were  found  to  be 
consistent  with  measurements  obtained  in  other  geophysical  flows. 

■  In  common  with  other  natural  flows  over  rough  boundaries,  the  bulk  of  Reynolds  stress  production  in 
Squaw  Creek  is  associated  with  infrequent  "events".  In  general  "events*  have  a  duration  greater  than  1.5 
seconds  and  a  return  period  greater  than  10  seconds. 

■  Typical  eddy  length  scales  determined  using  wavenumber  spectra  and  autocorrelation  functions  were 
found  to  be  0.52,  0J4  and  0.11m  for  u,  v  and  w  respectively. 

■  Spectral  slope  values  in  the  range  -1.85  to  -2^8  demonstrate  that  fully  isotropic  conditions  were  not 
detected  using  the  present  sensors.  This  is  considered  to  arise  through  a  combination  of  circu instances 
including:  selection  of  record  length;  digitization  rate;  low-pass  filler  cut-off;  and  the  physical  size  of  the 
EM  CM  heads.  The  individual  significance  of  any  single  factor  could  not  be  determined. 

■  Bed  shear  stress  values  obtained  from  EM  CM  measurements  using  TKE  methods  were  found  to  be 
consistently  lower  than  those  determined  using  RS  method.  Further,  location  of  the  EM  CM  sensors  well 
outside  the  "law  of  the  wall"  region,  high  bed  roughness  values  and  factors  outlined  above  probably 
contribute  to  the  disparity  between  u.  values  obtained  using  these  methods  and  some  of  those  obtained 
using  the  profile  method.  In  general,  however,  u.  values  obtained  here  are  in  good  agreement  and  are 
thought  to  be  reliable  indicators  of  bed  shear  stress  in  Squaw  Creek. 

■  The  drag  coefficients  and  roughness  length  values  obtained  using  the  present  measurements  of  bed 
shear  stress  and  flow  speed  are  found  to  be  consistent  within  data  series;  whether  calculated  from  the 
turbulence  data  series  or  from  the  velocity  profiles.  However,  bed  roughness  lengths  estimated  from 
turbulence  data  differ  from  those  estimated  from  profile  data. 

■  Roughness  length  values  obtained  from  near-bed  velocity  profile  data  scale  with  the  grain  size  of  the 
bed  material,  whilst  roughness  length  values  obtained  from  profile  data  which  includes  some  outer-layer 
observations  scale  with  the  dimensions  of  largest  elements  of  the  bed  armour. 

■  Detailed  velocity  profiles  taken  through  the  constant  stress  layer  and  the  outer-layer  indicate  the 
possible  presence  of  an  inflection  above  the  near-bed  layer.  The  presence  of  these  compound  velocity 
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profiles  may  indicate  the  presence  of  two  scales  of  bed  roughness  and  account  in  part  for  the 
discrepancies  noted  in  estimated  values  obtained  using  different  instrumentation  systems  and 
procedures. 

■  Manual  generation  of  SGN  at  0.5m  from  the  hydrophone  using  a  range  of  selected  pebble  sizes 
produced  acoustic  energy  significantly  higher  than  background  in  Squaw  Creek.  SGN  generated  in  this 
way  is  easily  detected  and,  based  upon  past  experience,  is  considered  to  be  typical  of  SGN  arising 
during  natural  bedload  transport 

■  The  spectra  for  the  bedload  detector  count  data  have  no  distinct  peaks.  Variance  in  contrast  is  spread 
fairly  evenly  across  all  the  frequencies  shown.  Consequently,  there  is  no  periodic  pulsing  during  the 
passage  of  the  sheets;  rather  transport  is  statistically  unsteady. 

■  The  statistical  distributions  of  width-summed  count  data  and  for  individual  sections,  at  both  the 
upstream  and  downstream  detectors,  conform  closely  to  positively-skewed  normal  distributions.  The 
positive  skew  may  originate  primarily  from  sampling  constraints  but  this  needs  further  investigation.  The 
similarity  between  sections  in  the  nature  of  the  distributions  and,  in  particular,  the  lack  of  bi-modality 
adds  weight  to  the  conclusion  that  there  is  no  pulsing  of  transport  and  no  movement  of  bedload  in 
distinct  and  spatially  discrete  bands. 

■  The  spectral  similarity  and  lack  of  phase  shift  between  upstream  and  downstream  series  and  between 
adjacent  sections  indicates  that  the  dynamic  response  of  the  bar-top  and  channel  talweg  are  in  all 
probability  contemporaneous  and  driven  by  a  forcing  function  which  operates  at  span-widths  equivalent 
to  the  channel  breadth  (10m)  and  the  reach  length  (30m). 

■  The  successful  application  of  a  general  statistical  transfer  function  to  sub-sets  of  the  upstream  and 
downstream  detector  count  data  indicates  that  coherent  sheet  structure  is  maintained  over  length  scales 
equivalent  to  the  experimental  reach  (30m).  In  other  words  the  sheet  retains  some  memory  of  upstream 
behaviour  as  it  translates  downstream. 

■  General  transfer  functions  may  be  suitable  to  model  sheet  dynamics  over  similar  length  scales  in 
other  gravel  bed  rivers  other  than  Squaw  Creek. 

■  The  finer  grain  sizes  from  the  mobile  sheet  art  trapped  in  the  interstices  of  the  bed  armour  as  the 
sheet  arrives  in  the  study  reach.  This  deposition  process  leads  to  physical  and  hydraulic  smoothing  of 
the  bed.  Towards  the  tail-end  of  the  sheets  passage,  the  finer  material  is  re-entrained  which  leads  to 
roughening  of  the  bed.  This  process  induces  a  response  in  hydraulic  parameters. 

■  In  any  modelling  exercise  the  effect  of  weak  bed  load  transport  on  bed  resistance  cannot  be  ignored, 
as  the  passage  of  discrete  bed  load  sheets  induces  strong  hysteretic  effects  in  hydraulic  parameters; 
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including  water  level.  The  inapplicability  of  stable-bed  equations,  such  as  Hey’s,  during  weak  bed  load 
transport  has  important  repercussions  for  the  use  of  general  flow  routing  models  for  shallow  mountain 
rivers. 

■  If  the  observations  and  conclusions  reported  here  also  apply  to  larger  gravel-bed  rivers,  then  there  are 
non-trivial  implications  for  calculations  required  for  engineering  management  For  example  maintaining 
adequate  freeboard  along  flood  defence  dykes,  positioning  of  water  abstraction  intakes,  and  navigation 
depths  depend  on  calculated  water  depths  being  accurate. 

4.1  Recommendations  for  further  research 

Field  Based  Research 


Additional  field  data  would  be  valuable  to  clarify  the  grain-sorting  mechanisms  as  bedload  sheets  pass 
and  the  interaction  with  the  bed.  In  particular  additional  data  at  higher  flow  levels  (than  are  reported 
here)  are  desirable.  These  data  can  be  obtained,  for  example,  using  a  vortex  sampler  or  Helley-Smith 
type  sampler,  together  with  fine-resolution  profiling  of  the  bed  elevations.  These  grain-sorting  data  need 
relating  to  hydraulic  parameters  calculated  born  detailed  near-bed  velocity  profiles  obtained  using 
electromagnetic  current  meters  as  well  as  near-bed  turbulence  data. 

The  magnetic  detector  system  and  the  acoustic  detection  system  provide  an  opportunity  to  obtain  long- 
runs  of  data  on  bedload  transport  at  a  temporal  resolution  necessary  to  relate  particle  dynamics  to  the 
scale  of  turbulent  fluctuations.  Not  only  does  this  have  an  important  research  context,  but  also  offers  the 
possibility  of  routine  monitoring  of  bedload  in  remote  locations. 

The  installations  at  Squaw  Creek  are  unique  and  provide  an  opportunity  to  obtain  the  integrated  data¬ 
sets  required  to  elucidate  the  mechanisms  of  coarse  bedload  sheet  dynamics  and  to  develop  further 
automated  bedload  monitoring.  It  is  recommended  that  consideration  should  be  given  to  developing  a 
further  research  and  technical  development  programme  at  Squaw  Creek. 

Modelling  Unsteady  Bed  load  Transport 

The  particle  interactions  with  the  flow  and  the  bed  material  can  be  modelled  using  parallel  processing 
techniques*4.  The  existing  detailed  comprehensive  hydraulic,  channel  geometry  and  bedload  data  could 
be  used  to  construct  a  bedload  transport  model  mimicking  unsteady  behaviour  based  upon  fundamental 
physical  principles.  The  implications  for  future  generations  of  engineering  models  such  as  HEC-2  should 
be  a  key  focus  of  the  investigation.  It  is  recommended  that  a  computer-based  mathematical  model 
should  be  developed  which  simulates  the  passage  of  coarse  bedload  in  gravel  rivers  as  unsteady  waves 
of  material 
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ANNEX  1:  TURBULENCE  DATA  PREPROCESSING  AND  ANALYSIS. 

1.0  Scope  of  the  work 

In  this  section  the  methodology  used  to  analyze  the  EM CM  data  is  detailed  and  problems  likely  to  affect 
data  quality  and  inteipretation  ate  highlighted.  In  general,  data  processing  aims  to  quantify  flow 
turbulence  and  Reynolds  stresses  and  to  relate  these  properties  to  bedload  dynamics. 

2.0  Data  Screening 

All  data  from  the  data  logger  are  downloaded  via  a  PC  to  a  work-station.  Single  files  containing  EMCM 
and  SGN  data  are  prepared  for  each  measurement  run  and  channel  numbers  are  assigned  for  each  data 
set.  Raw,  uncalibrated  data  are  then  plotted  and,  in  order  to  assess  data  quality,  the  mean  and  standard 
deviation  for  each  data  channel  are  calculated. 

3.0  Data  editing 

In  order  to  screen  and  edit  data  sets,  files  are  processed  using  interactive  software,  (SERPL04).  This 
procedure  involves  the  identification  of  infirequert  data  spikes  by  eye  and  subsequent  ’flagging"  of  all 
suspect  points  in  a  given  time  series.  These  flags  are  easily  picked  up  in  later  processing  and  are 
removed  using  linear  interpolations  between  adjacent  data  points.  Use  of  SERPL04  also  facilitates  first 
look  inter-comparisons  between  data  channels  and  a  rapid  "zoom-in”,  "zoom-out"  and  scan  option 
enables  rapid  assessment  of  data  quality. 

441  Offsets  and  trend  removal 

Zero  offset  errors  were  present  on  all  channels.  These  were  measured  carefully  in  the  field  before  and 
after  individual  test  runs.  The  resulting  flume  calibrations  and  offsets  are  applied  to  the  appropriate  data 
channels  in  each  burst  record. 

A  common  problem  with  output  from  all  EMCM  systems  is  a  tendency  for  the  zero  level  to  drift 
through  time.  However,  in  the  case  of  data  described  here,  the  zero  drift  was  very  low  during  test  runs 
and  could  be  removed  effectively  using  a  least  squares  trend  removal  technique. 
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5.0  EMCM  axis  rotation 

When  deployed,  great  care  was  taken  to  ensure  the  EMCM  axis  was  in  line  with  the  principle  streamline 
of  Squaw  Creek.  However,  as  accuracy  of  this  technique  is  subjective,  it  was  necessary  to  calculate 
EMCM  head  rotation  angles  for  perfect  alignment  with  the  flow.  This  is  achieved  by  plotting  the 
apparent  vertical  mean  velocity  versus  the  mean  horizontal  velocity  for  successive  blocks  in  a  given 
record.  Each  plot  shows  a  general  trend  corresponding  to  the  mean  tilt  angle  throughout  the  deployment, 
(correlation  coefficients  >  95%  significant).  Uncertainties  in  zero  drift  removal,  variations  in  sensor 
attitude  relative  to  the  flow  and  bed  through  time  and  calibration  errors  combine  to  give  rise  to  the 
scatter  of  data  values  in  such  plots.  Following  the  present  configuration  of  the  EMCM  beads,  application 
of  a  horizontal  rotation  through  45„  ($)  is  also  required  for  the  channels  containing  the  U  and  V  flow 
components.  This  procedure  gives  time  series  of  the  turbulent  flow  components  aligned  along  the 
fore/aft  and  transverse  axes  of  the  current  meters.  Rotation  procedures  ate  described  in  6.3. 

63  Data  Analysis 

6.3.1  Summary  of  EMCM  data 

In  the  processing  of  EMCM  data,  both  time  series  and  summary  statistics  are  generated  in  a  form  which 
facilitates  direct  comparison  between  experiments.  In  the  description  of  calculated  quantities  given  below 
the  following  notation  is  used: 

z  Height  of  instrument  above  the  bed,  (no  correction  for  local  scour  applied). 

U  Mean  velocity  along  the  fore/aft  EMCM  axis,  (positive  for  downstream  flow). 

v  Mean  velocity  along  the  transverse  EMCM  axis,  (positive  for  flow  from  port  to 

starboard  ie;  right  bank  to  left  bank  looking  downstream). 

u  Horizontal  turbulent  velocity  component. 

v  Transverse  •  •  • 

w  Vertical  "  '  ' 

p  Fluid  density 
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6.3.2  Spectra 

Eneigy  spectra  of  u,  v  and  w  arc  obtained  from  the  9.1  minute  runs  (16  384  data  values)  using  the  fast 
Fourier  transform  (FFT)  routine.  The  raw  spectral  estimates  arc  smoothed  using  a  weighted  avenging 
scheme  and  plotted  against  frequency  on  log  axes  (  a  form  generally  used  for  turbulence  spectra).  The 
cospectn  uw  and  vw  arc  obtained  by  combining  raw  spectral  estimates  and  subsequently  smoothed 
using  a  weighted  scheme.  Cospectra  arc  plotted  in  the  conventional  "equal  area-equal  energy  form". 

€33  Time  series  and  statistics 

As  described  above,  EMCM  data  requires  suitable  rotation  in  both  the  horizontal  and  vertical  planes. 
These  operations  are  defined  below  where  the  subscripts  p  and  s  indicate  port  and  starboard  EMCM 
heads  respectively. 

vertical  rotation  angle,  0,  =  tan1(w(/uI1)  ; 

*  "  *  02  =  tan'l(w,/u1)  ; 

applying  vertical  rotation,  wp  =  wpcos  Oj-u^sin  0, ; 

■  "  *  u,,  =  w^in  0,-UpCos  0,  ; 

*  "  "  w,  =  waco6  02-ujsin  02 ; 

"  "  •  u,  =  w^in  02-u,cos  02 ; 

applying  horizontal  rotation,  U  =  (ujsin  4>)  +  (ubcos  d>)  ; 

"  "  "  V  =  (u/in  d>)  -  (ubcos  <l>)  : 

Following  rotation,  software  analyzing  EMCM  data  produces  useful  summary  statistics  (see  examples  in 
Table  3).  In  addition,  the  following  30  Hz  time  series  are  calculated. 


current  speed, 


S  =  (U*+V?)w  ; 
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current  direction 


lPe=  tan*1  (V/U),  degrees  anti-clockwise  of  the  foie/aft  axis  ; 


fore/aft  Reynolds  stresses, 


Jtej  *  ~vw. 


transverse  Reynolds  stresses. 


*  ~VW, 


stress  direction. 


t,  -  tan"  ((-v*y(-«w)); 


stress  direction  re.  current. 


*  V*c  ; 


stress  magnitude, 


T/p  *  ((-iov*)+(-vw*)w; 


and,  turbulent  kinetic  energy, 


E  *  1/2 (Zf+vf+w?). 


63A  Errors 

Whilst  every  precaution  is  made  to  exclude  errors  from  the  processed  data  some  unwanted  "noise"  must 
be  tolerated.  Error  in  EM  CM  rotation  angle  give  rise  to  large  errors  in  Reynolds  stress  values.  Using  the 
criteria  outlined  by  Souls  by70,  it  is  considered  that  although  suffering  from  the  same  minor 
imperfections,  the  processed  data  are  acceptably  accurate  and  ate  representative  of  turbulent  flow 
conditions  in  Squaw  Creek. 


Table  1 


Example  of  particle  size  distribution  of  bed  material  (mm). 
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Table  2 


Summary  of  U,  RMS  values  for  u,  v  and  w,  u.  values  from  TKE  and  RS  measurements  and  Cd  and  zD 
values  during  runs  4  to  8  in  Squaw  Creek,  June  1992. 
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Table  3 


Useful  data  summary  statistics  output  for  Run  4,  Squaw  Creek  data. 
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Figure  5. 


Figure  6. 


Schematic  diagram  of  1981  detector  (after  Custer  &  Ergenzinger,  1983) 


Schematic  diagram  of  the  1986  detector  (Irom  Spieker  &  Ergenzinger,  1990) 
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Figure  7. 
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Figure  13. 
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Figure  14. 
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Figure  18. 
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Figure  24. 
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Figure  26. 
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Figure  27. 
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Figure  29. 
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